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We present numerical results on the capacities of two-qubit unitary operations for creating entan- 
glement and increasing the Holevo information of an ensemble. In all cases tested, the maximum 
values calculated for the capacities based on the Holevo information are close to the capacities based 
on the entanglement. This indicates that the capacities based on the Holevo information, which 
are very difficult to calculate, may be estimated from the capacities based upon the entanglement, 
which are relatively straightforward to calculate. 
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I. INTRODUCTION 

A nonlocal operation is one which operates on distinct 
subsystems, and can not be expressed as a tensor prod- 
uct of operations on the individual subsystems. Such 
operations arise via an interaction Hamiltonian between 
the subsystems. Nonlocal operations may be used to cre- 
ate entanglement between subsystems, and also to per- 
form classical communication. In fact, it is not possible 
to achieve either of these tasks without an interaction 
between the subsystems. In characterizing quantum op- 
erations, it is thus important to determine the capac- 
ities for creating entanglement or performing communi- 
cation. These capacities provide a characterization of the 
strength of the operations [l| • 

Shared classical information may be considered to be 
the classical equivalent of entanglement. Therefore it is 
reasonable to consider the process of classical commu- 
nication to be the classical equivalent of entanglement 
creation. One may therefore expect that there is a close 
relationship between the capacities for these two tasks. 
In fact, it has been shown i^^^l] that there are a number 
of inequalities between the various capacities for classi- 
cal communication and entanglement creation. In this 
paper we make a direct numerical comparison between 
these capacities. 

The capacities of operations for creating entanglement 
have been studied extensively H S 0, EM E 111 Hi- 
lt is relatively straightforward to determine the entan- 
glement capability for infinitesimal operations 0- For 
finite operations, most results are restricted to numerical 
results for certain classes of two-qubit operations [lol |. 
We study the same classes of operations here, and the 
results we present for the entanglement reproduce those 
given in Ref. T^, except for some data points where we 
believe our results arc more accurate. 

The capacities for classical communication were ini- 
tially considered for simple operations, such as the CNOT 
and SWAP operations 13, 14] . It is more difficult to con- 
sider the classical communication for general operations, 
because general operations can not be used to perform 



perfect communication. Ref. '2'| introduced asymptotic 
capacities, where the average communication when the 
operation is performed a very large number of times is 
considered. When the operation is performed a large 
number of times, it is possible to use error correcting 
techniques to reduce the probability of error to be arbi- 
trarily small. 

For the purposes of numerical investigation, it is not 
practical to use asymptotic capacities, as the number of 
variables that would need to be optimized over to ob- 
tain a reasonable approximation is extremely large. A 
far more practical type of capacity is that based on the 
Holevo information of ensembles. For these capacities, 
only the states and probabilities in the initial ensemble 
need be optimized over in the numerical analysis. In 
addition, there are connections between the capacities 
defined via the change in Holevo information and the 
asymptotic capacities 2]. 

Although it is feasible to calculate capacities based on 
the Holevo information, it is far more computationally 
difficult than evaluating the entanglement capacities. It 
would therefore be useful if it were possible to use the 
entanglement capacities to estimate the capacities based 
on the Holevo information. In this paper we numerically 
study the relationship between these capacities, and show 
that although they are not equal, the values we have 
calculated are quite close. 

In Sec. ^ we explain the definitions of the various ca- 
pacities. In Sec. mil we give numerical results for the 
communication capacity based on the Holevo informa- 
tion obtained when the initial ensemble has zero Holevo 
information, and compare these capacities to the entan- 
glement that may be created from initial states that have 
zero entanglement. Then in Sec. IIVI we give results for 
the increase in Holevo information for general initial en- 
sembles, and compare these capacities to the increase in 
entanglement for arbitrary initially entangled states. We 
conclude in Sec. and in the Appendix give a detailed 
explanation of the numerical techniques used to calculate 
the results presented. 
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II. DEFINITIONS 

First we summarize the definitions of the various ca- 
pacities. Throughout this paper we divide the system 
into two subsystems, A and B, and denote the Hilbert 
spaces by Ha and Hb- The party in possession of sub- 
system A will be referred to as Alice and the party in 
possession of subsystem B will be referred to as Bob. 
The subsystems A and B are divided into further sub- 
systems: 

Ha - Hau ® Wa.„, , Hb = Ubu ® Ws.„c • (1) 

The operation U acts only upon "Hau ® Hbui and the 
Hilbert spaces Ti-A^na and "Hb^^,, are ancillas. The ancillas 
may have dimension that is arbitrarily large but finite. 

There are two main ways of defining capacities for en- 
tanglement. The first is the entanglement that may be 
obtained when the initial state is unentangled 

Eu= sup E{U\4>)a\x)b). (2) 

\4>)AenA,\x)BenB 

The quantity E{- ■ ■) is the entropy of entanglement 
Em) = SlTiAimmi where Sip) = -Tr(plogp). 
Throughout we employ logarithms to base 2, so the en- 
tanglement is expressed in units of ebits. The second 
definition is the maximum increase in entanglement when 
the initial state may be an arbitrary pure entangled state. 

AEu= sup (3) 

It is also possible to define the capacity in terms of the 
entanglement of formation and allow mixed states; we do 
not separately consider this case, because allowing mixed 
states does not alter the capacity 

Another alternative definition of the entanglement ca- 
pacity is based on the average entanglement that may be 
obtained in the limit that the operation is performed an 
extremely large number of times and the initial state is 
unentangled |3- It is shown in Ref. that this capac- 
ity is equal to the maximum increase in entanglement as 
defined in Eq. (jSJ. Thus the numerical results presented 
also apply to the asymptotic entanglement capacity. 

In the numerical search, it is not possible to consider 
ancilla spaces with arbitrarily large dimension. It is 
known jl| that the ancilla spaces need have dimension 
no larger than the Hilbert spaces Hau and Hbu fo^' the 
capacity with initially unentangled states. It has also 
been found numerically 10] that the same is true for the 
capacity where initially entangled states are allowed. 

In the results presented below, we use a fixed dimen- 
sion for the ancilla spaces. For simplicity we use equal 
dimensions on the two ancilla spaces Ha^^^ and 7is^„^. 
We use a superscript on the capacity when it is neces- 
sary to indicate the dimension of the ancilla space used. 
For example, AiJ^f"* is the maximum change in entangle- 
ment when the ancilla spaces are each of dimension 4. 



When we refer to multiple results with different ancilla 
dimensions we will use a superscript (*). We will omit 
the superscript in the case of Eu when the dimensions of 
the ancilla spaces are at least as large as those of Hau 
and Ti. Bu , because this is known to be sufficient to obtain 
the capacity for arbitrarily large ancilla. 

The classical communication capacities that we will 
consider are based upon the Holevo information of en- 
sembles. An ensemble is a set of states {|$i)A_B} that 
are supplied with probabilities pi. Each state \^i)AB 
is a pure state shared between Alice and Bob, and Al- 
ice chooses the index i. The ensemble is denoted by 
£ — {pi,\^i)AB}- We also define the ensemble of reduced 
density matrices possessed by Bob as 

E = TTA£^{p^,P^}, (4) 

where pi = Ti-A\^i)AB{^i\- The Holcvo information of 
the ensemble E is given by 

X{E) ^ S (^p^^p}j -J^P^Sip.). (5) 

From the Holevo-Schumacher- Westmoreland theorem 
ITgI I , the Holevo information gives the average com- 
munication that may be performed from Alice to Bob by 
coding over multiple states. 

Similarly to the case for entanglement, we may define 
capacities based on the maximum change in Holevo in- 
formation. One definition that we will use is the maxi- 
mum final Holevo information when the initial ensemble 
has zero Holevo information. For the initial ensemble, 
we have an initial state \'4')ab, and Alice encodes i by 
applying a local operation V^^K For the capacity, the 
supremum is taken over the initial state |V')as, the en- 
coding operations V'-*' and the probabilities p^: 

XU= sup x(p^,^rAUV^'^iJ)AB), (6) 

p.,y(').|V>AB ^ ^ 

where Tr^|0) = Tia\4'){4'\- Note that the notation we 
are using here differs from Ref. 0, where the symbol 

^Xc/'"^ was used for this capacity. 

We may also define the maximum change in Holevo 
information when the initial ensemble is arbitrary: 

Axu = supxi^^AUS) - x(Tr^£). (7) 

£ 

Here we are using the notation conventions 

C/£ = fe,C/|$,)}, (8) 
Tr^£ = fe,Trx(|$,))}. (9) 

This capacity is equivalent to the capacity A^^^'*-* de- 
fined in Ref. As shown in Ref. 2J, this capacity is 
equal to the average entanglement-assisted communica- 
tion that may be performed from Alice to Bob. Therefore 
this quantity may be interpreted as the asymptotic com- 
munication capacity, just as AEu may be interpreted as 
the asymptotic entanglement capability. 
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One may also interpret xu in terms of asymptotic ca- 
pacities. The capacity xu gives the Holevo information 
after a single application of the operation U. This com- 
munication can not actually be performed for a single 
ensemble; it is necessary to code over multiple states to 
perform this average communication. Therefore xu niay 
be interpreted as the asymptotic communication capacity 
if the communication protocol is limited to the relatively 
simple scheme where coding is performed over multiple 
final states. This is equivalent to restricting all the appli- 
cations of U to be performed at the same time (on input 
states that are not entangled with each other), rather 
than allowing the output of one application of U to be 
used as part of the input to another application of U, as 
in the general case. 

As in the case of the entanglement, it is not possi- 
ble to use ancilla spaces of arbitrarily large dimension in 
the numerical search. In addition, it is not possible to 
use arbitrarily large numbers of states in the ensemble. 
In the results presented below, we perform calculations 
for restricted ensembles where there is a fixed number 
of states in the ensemble and a fixed dimension for the 
ancilla spaces (the ancilla spaces are again taken to be of 
equal dimension). We use superscripts on the capacities 
to indicate the number of states in the ensemble and the 
dimension of the ancilla spaces. For example, Axl^'^"* 
is the maximum change in Holevo information for two 
states in the ensemble and ancillas each of dimension 4. 
We use a superscript (*) to refer to multiple capacities 
with different ancilla dimensions or ensemble sizes. It 
must be emphasized that our use of superscripts in this 
paper differs from that in Ref. 0- 



III. CAPACITIES FOR ZERO INITIAL HOLEVO 
INFORMATION 

It is clear that the capacity xu is an analogous quantity 
for communication to Eu for entanglement. Similarly 
Axu is analogous to AEu for entanglement. In this sec- 
tion we perform a direct numerical comparison between 
the two capacities xu and Eu. In the next section we 
compare the capacities Axu and AEu. 

In this paper we concentrate on two-qubit operations. 
It is not possible to perform calculations for the entire 
range of two-qubit operations. To make the problem fea- 
sible, we only consider a limited number of examples of 
two-qubit operations. In particular, we consider opera- 
tions of the form: 

C/i(a) = [/d(a,0,0), (10) 
U2{a) ^Ud{a,a,0), (11) 
U3{a) = Ud{a,a,a), (12) 

where 



Here CTj are the Pauli sigma operators. The operations 
Ui, U2 and C/3 correspond to the CNOT, DCNOT and 
SWAP families of operations considered in Ref. . 

In order to consider the complete range of two-qubit 
operations in the case of entanglement, it would only be 
necessary to consider operations of the form H13() . with 
7r/4 > ai > a2 > 03 > |^. This derivation relies on 
the fact that any two-qubit operation may be simplified 
to one of the form (|13l l with 7r/4 > ai > ±a2 > > 
using local operations 0, 0, 0, In addition to using 
local operations, the derivation in Ref. |3| relies on the 
fact that the entanglement capabilities of U and U* are 
identical (which implies that all the ai may be taken to 
be positive). 

Similarly, for the Holevo information, x(TrA^) = 
x(TrAf *) and xC^^aUS) = x(TrA [/*£•*). Thus the ca- 
pacities of U and U* to increase the Holevo information 
are identical. Therefore, in order to obtain results for the 
complete range of two-qubit operations, it would only be 
necessary to consider operations of the form H13|l with 
7r/4 > Q!i > a2 > CX3 > for both entanglement and 
Holevo information. This restriction on the values of the 
ai defines a three dimensional region of values. In this 
paper we do not consider the entire region; however, the 
operations Ui, U2 and C/3 which we consider form three 
lines on the boundaries of this region. 

Sophisticated numerical maximization techniques were 
used to find the values of Eu and Xu^\ the details are 
given in the Appendix. The numerical results for Eu-^ 
and Xu^ are shown in Fig. ^ The values of Eu-^ were 
determined in intervals of 7r/400 for a, and are shown 
by the solid line. An ancilla dimension of 2 is sufficient 
to obtain the asymptotic capacity Eu for all two-qubit 
operations In all the calculations presented for Eu 
in this section, an ancilla dimension of 2 was used; we 
have therefore omitted the superscript on this capacity. 
In addition, it was found that, for Ui, the same capacity 
was obtained without ancilla, in agreement with Ref. |9| . 
(*) 

The values of Xjj^ were determined in intervals of 7r/40 
for a, with ancilla dimensions of 1, 2 and 4, and ensembles 
with 2 or 4 members. An ancilla dimension of 1 is equiva- 
lent to no ancilla; we will take the case of no ancilla to be 
an ancilla dimension of 1 for convenience in the notation. 
Just as in the case of the entanglement, these results 
indicate that the maximum final Holevo information is 
obtained without ancillas. In addition, only two states 
in the ensemble are required. We find no improvement 
in using up to four states in the ensemble, and ancillas 
of dimension up to four. It was found that the max- 
imum increase in the Holevo information was obtained 
with equal probabilities when there were two states in 
the ensemble. When four states in the ensemble were 
used, in most cases two of the probabilities approached 
zero whereas the other two approached 1/2. 

In addition, note that there is no difference between the 
results obtained for the maximum final entanglement and 
the maximum final Holevo information. It was found that 
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FIG. 1: The increase in entanglement for no initial entangle- 
ment, Eui, and the increase in Holevo information for zero 
initial Holevo information, Xjji ■ The values of Eu^ are shown 
as the solid line, and the values of Xul^^ ' Xuf^' Xuf^ ^^'^ 
Xuf^^ are shown as the circles, crosses, plusses and squares, 
respectively. These symbols overlap and are not visible sepa- 
rately. 

the two capacities agreed to precision better than one 
part in 10^'*. Our results do not prove that Eu-^ — xun 
because it is not possible to test ensembles vifith arbitrar- 
ily large numbers of states or allow the ancilla dimension 
to be arbitrarily large. Nevertheless, our results strongly 
indicate that Ejji is equal to xui ■ 

Numerical results for Eu^ and Xu2 shown in Fig. |21 
As for Ejj-^ , the value of Ejj^ was determined in intervals 
of 7r/400 for a, and an ancilla dimension of 2 was used in 
order to obtain the asymptotic entanglement capability 
Ejj^ . The values of ^^^^ determined in intervals of 
7r/40 for a, and with ancilla dimensions of 2, 3, 4 and 5. 
In each case the ensemble consists of four states. 

In contrast to the case for the entanglement, there are 
improvements in using ancilla dimensions higher than 2 
when considering the Holevo information. For an ancilla 
dimension of 2, the value of Xc/j^"* is equal to Ejj^ for the 
lower values of a, but is less than Eu^ for the upper values 

of a. For an ancilla dimension of 3, the values of Xj/^'^'' 
are larger, but still less than Ejj^ towards the right side 
of the plot. For an ancilla dimension of 4, Xuf^^ equal 
to Ejj^. When the ancilla dimension is increased to 5, 

the capacity Xuf'^ is no larger than Xjjf^^- These results 
indicate that an ancilla dimension of 4 may be sufficient 
to achieve the capacity XU2 j and that this capacity is 
equal to Ejj^ ■ 

Similar calculations were performed with equal prob- 
abilities for each of the states in the ensemble. In most 
cases the same results were obtained as when arbitrary 
probabilities were allowed. However, for one data point 
for Xuf^ (indicated by the circle in Fig. |2Jl, a larger 
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FIG. 2: The increase in entanglement for no initial entangle- 
ment, Eu2, and the increase in Holevo information for zero 
initial Holevo information, Xu2 ■ values of Eu2 are shown 

as the solid line, and the values of Xu2^^ y Xc/j'^' and Xc/j''' are 
shown as the plusses, crosses and squares, respectively. All 
results for are equal to Xc/j'''') and are not shown be 

separate symbols. The data point indicated by the circle is 
for xl/j^"* when the probabilities are restricted to be equal. 

Holevo information was obtained when unequal proba- 
bilities were allowed. The numerically derived optimal 
ensemble had three members with non-negligible prob- 
abilities: two probabilities were approximately 41.39%, 
and one was 17.22%. 

Results for Eu^ and Xu^ ai'e shown in Fig. |31 Again 
the value of Ejj^ was determined in intervals of 7r/400 for 
a, and with an ancilla dimension of 2. The values of xus 
were determined for ancilla dimensions of 2, 3, 4 and 5, 
and with four states in the ensemble. 

As in the case of Xu2^\ the values of Xuf' are signifi- 
cantly below Ejj^ for most of the plot. When the ancilla 

dimension is increased to 3, the values of Xc/J'^'' are no- 
ticeably larger, and closer to Ejj^. In addition, three of 
the values, at a = tt/40, 27r/40 and 37r/40, are actually 
larger than E'f/g (see the subplot in Fig. I^)). The difference 
is small, less than 0.02, but it is sufficient to demonstrate 
that Ejj^ is not equal to xua ■ For an ancilla dimension 
of 4, Xua'^^ i^ ^*^t smaller than Eu^ for any of the data 

points. In addition, Xuf^^ larger than Ejj^ for the same 
three values of a as for an ancilla dimension of 3. In fact, 
Xuf^^ i^ equal to Xuf ^'^^ these three data points. When 
the ancilla dimension is increased to 5, Xj/J^'' is no larger 

than Xu'f^ ' again indicating that an ancilla dimension of 
4 is sufficient to achieve the capacity xua ■ 

To summarize, we have found that, for the operation 
Ui, ancillas do not increase the capacity xu, just as in 
the case for the entanglement. For the cases U2 and U3, 
the capacity xu increases with ancilla dimension up to an 
ancilla dimension of 4, but is unchanged when the ancilla 
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FIG. 3: The increase in entanglement for no initial entangle- 
ment, Eu^, and the increase in Holevo information for zero 
initial Holevo information, Xua- "^^^ values of Eu^ are shown 

as the solid line, and the values of X[/3^\ ^c/a'^' '^t^J*' 
shown as the plusses, crosses and squares, respectively. All 

values of identical to Xu '^\ S'lid are not shown in- 

(*) 

dependently. The subplot shows the difference Xu^ ~ Eu^ for 
Xuf^ (crosses) and Xu^^^ (squares). 

dimension is further increased to 5. This is in contrast 
to the case for the entanglement, where the maximum 
entanglement capacity Eu is obtained for an ancilla di- 
mension of 2 13] ■ 

The results indicate that xu is equal to Eu for the 
operations Ui and U2, though it is possible that Xjj is 
increased for larger ancilla dimensions or ensemble sizes. 
For the operation U3, it is possible to obtain slightly 
higher values of Xu\ demonstrating that xu is not equal 
to Eu for this operation. In all cases tested, we find that 
X(7 > Eij. In only a small number of cases have we found 
that Xu 7^ Eu, and in these cases the differences found 
are only small, suggesting that Eu is an excellent, and 
efficient, estimator oi xu- 



IV. CAPACITIES FOR ARBITRARY INITIAL 
ENSEMBLES 

Next we consider the capacities Axu and AEu- These 
capacities are more general, in that arbitrary initial states 
or ensembles are allowed. Analytic results for the relation 
between these capacities have been derived in Ref. . In 
this reference it is shown that, for two-qubit operations, 
Axu > AEu ■ If AEu may be achieved with a particular 
ancilla dimension, then an increase in Holevo information 
equal to AEu may be achieved with the same ancilla 
dimension, and with four members in the ensemble. 

In principle it is possible that there is no ancilla dimen- 
sion that achieves AEu , and instead AEu is approached 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 



a 

FIG. 4: The increase in entanglement for arbitrary initial 
state, AiS^f*' , and the increase in Holevo information for arbi- 
trary initial ensemble, Axij^- The values of /^E^^ and AE^^ 
are shown as the solid line, and the values of Axij^^^ and 
Ax^f^ are shown as the circles and plusses, respectively. 

in the limit of large ancilla dimension. However, in prac- 
tice it has been found that, for two-qubit operations, it 
appears to be possible to achieve AEu with an ancilla 
dimension of 2 This means that it should be possi- 
ble to achieve an increase in Holevo information equal to 
AEu with an ancilla dimension of 2. 

The capacities Axu^ and AE^J^ are shown in Fig. ^ 

Each capacity was determined in steps of 7r/40 in a. It 

(*) 

was found that the entanglement capacity AE\j^ did not 
increase beyond that for no ancilla as the ancilla dimen- 
sion was increased up to 5, in agreement with the result 
given in Ref. Throughout this section we use the 

superscript asterisk when referring to results for AEfj , 
because it has not been proven that the asymptotic ca- 
pacity is achieved for an ancilla dimension of 2. 

The capacities Ax[jl without ancillas and with ancillas 
of dimension 2 are shown in Fig. 0] In both cases these 
capacities are for ensembles with two states. It is found 
that, even without ancilla, the capacity Axul greater 

than the values calculated for Ai?^^ . The only cases 
where there is equality are the trivial cases where a = 
or 7r/4. These results indicate that there are operations 
for which there is the strict inequality Axu > AEu- 

In addition, the capacity Axu^ is slightly increased by 
adding an ancilla. This is not so easily seen in Fig. m to 
make this difference visible, the differences between the 
capacities Axu^ with ancilla and the capacity with no 
ancilla Axu^ are plotted in Fig. |31 It can be seen that 

there is a small but significant increase in A^^^ when an 
ancilla is allowed. 

In addition, there is a further improvement in using an 
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FIG. 5: The increase in Holevo information for arbitrary ini- 
tial ensemble, Ax^^, relative to that for the case where the 
ensemble has two members and there is no ancilla, Ax^'^'^' 
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FIG. 6: The increase in Holevo information for arbitrary ini- 
tial ensemble, Ax^^, relative to that for the case where the 
ensemble has two members and the ancilla is of dimension 2, 
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The values of Axjjf^ 
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Ax^/i - Ax^; 



^jj^ — "A(7^ Bjce shown as the crosses, squares and 
asterisks, respectively. 



ancilla dimension of 3 rather than an ancilla dimension 
of 2. There are further increases in the capacity as the 
ancilla dimension is increased to 4 and 5, as shown in 
Fig. IHl These results indicate that the true asymptotic 
capacity Axui is not actually achieved for any particular 
ancilla dimension. However, each increase in the capacity 
with the ancilla dimension is smaller than the previous, 
indicating that the results calculated here should give a 
good approximation of Axui ■ 

Calculations were also performed for ensembles with 4 
members, and ancilla dimensions up to 4. It was found 
that, for all ancilla dimensions tested, there was no in- 
crease in the capacity when the number of states in the 
ancilla was increased. In addition, it was found that there 
was no increase in the capacity with up to 8 states in the 
ensemble and no ancilla. 

The results for Ax[;*^ and AeIJ^ are shown in Fig. [3 
In each case shown, ensembles with four states were used. 
In the case without ancilla, it was found that Ax^^J^-* and 
AE^j^ were equal. When an ancilla is included, there 

is a significant increase in both Ax^^' and AE'J^ . In 
particular, these have a maximum of 2, rather than 1 as 
in the case without ancilla. 

Note also that the value of AE^^ is increased when 
an ancilla is added for each of the values of a except the 
trivial points at a = and 7r/4. In contrast, for the 
data shown in Ref. |0 there was no visible increase in 

AE)j^ when the ancilla was included for another three 
data points (at a = 7r/40, 27r/40 and 37r/40). We sus- 
pect that this is because the optimization found the local 
maximum corresponding to the solution with no ancilla, 
rather than the global maximum. 
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FIG. 7: The increase in entanglement for arbitrary initial 
state, AE\^^ and the increase in Holevo information for arbi- 
trary initial ensemble, Ax^^ ■ '^^^ values of AE^^ and AE^^ 
are shown as the dotted and solid lines, respectively, and the 
values of Ax^^'^' , AXf/^^^ and Axjj^^^ are shown as the circles, 
plusses and crosses, respectively. 

We found that using ancilla dimensions above 2 up 
to an ancilla dimension of 5 did not increase AEjj^ , in 
agreement with Ref. id . When the ancilla was included, 

Axu^ was slightly greater than AEjj^, just as in the case 

of the operation Ui. In addition, it was found that Ax\j^ 
was further increased as the ancilla dimension was in- 
creased beyond 2 (see Fig. (S)). In this case the difference 
is somewhat greater, being around 0.02 rather than 10"*^, 
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FIG. 8: The increase in Holevo information for arbitrary ini- 
tial ensemble, A^^^^, relative to that for the case where the 
ensemble has four members and the ancilla is of dimension 2, 

Ax^^'^'. The values^of Ax^'f " ^x'uT^ ^x'uT " ^x'uf 
and Ax^f' — Ax^^^"* are shown as the crosses, squares and 
asterisks, respectively. 




FIG. 9: The increase in entanglement for arbitrary initial 
state, AEurf , and the increase in Holevo information for arbi- 
trary initial ensemble, Axua- The values of AEu^ are shown 
as the solid line, and the values of Ax^^'^' , Ax^^j^' and Ax^^'*' 
are shown as the plusses, crosses and squares, respectively. 
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but the values still appear to be converging for large an- 
cilla dimension. It is expected that the results shown are 
a good approximation of the true value of Axu2 ■ 

The results for the third operation, U^, are shown in 
Fig.|51 All results here are for ensembles with four states. 
In this case it was found that, if the ancilla had dimen- 



sion 2, the values of Axfjf'^ and AE)J^ were identical. 
In other respects the results were similar to those for the 
operation U2- The value of i^E\j^ was not increased by 
increasing the ancilla dimension above 2, as for the op- 
erations Ui and C/2. The value of ^Xui, increased 
for larger ancilla dimensions, so for these larger ancilla 
dimensions ^^Xui equal to i^E\j^ . 

The difference between ^Xui, higher ancilla dimen- 
sions and Axj^J^'' is shown in Fig. ^| In this case the 



p(2) 



difference between AxJ^'''^ 



W3 



and l^Xut^ 



is almost 0.05, 



which is larger than for both Ui and U2- Nevertheless, 

the maximum difference between l^x^uf^ ^'^'^ ^-^ul 
still comparable with the results for Ui and C/2. 
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FIG. 10: The increase in Holevo information for arbitrary 
initial ensemble, Ax^^'' , relative to that for the case where the 
ensemble has four members and the ancilla is of dimension 2, 



Ax[,'''' . The values of Ax[,'f' - 



-Ax'^f^nAAx'^f^ 



are shown as the crosses and squares, respectively. 



CONCLUSIONS 



To summarize, our results strongly indicate that Axu 
is strictly greater than AEu for most two-qubit opera- 
tions, rather than simply greater than or equal to, we 
stated in Refs. |^, In addition, our results show that 
the maximum change in Holevo information is not ob- 
tained for ancilla dimensions of two, as appears to be the 
case for the entanglement. The indications are that there 
is no ancilla dimension that is sufhciently large that Axu 
is achieved, so Axu would only be obtained asymptoti- 
cally in the limit of large ancilla dimension. 



Our results demonstrate that, for two-qubit opera- 
tions, there are close relationships between the capaci- 
ties based on the Holevo information and the capacities 
based on the entanglement. In most cases, the largest 
values calculated for xu (with the largest ancilla dimen- 
sion and ensemble size) were equal to Ejj . In some cases 
values of xu were calculated that are slightly above Ejj, 
but the difference was very small. In addition, the maxi- 
mum values calculated for Axu were very close to those 
for AEjj. These results indicate that the communication 
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capacities xu a-nd A^tj, which are very difficult to calcu- 
late, may be estimated from the entanglement capacities 
Eu and AEu, respectively, which may be calculated far 
more rapidly. 

We have also found that, in all cases tested (except 
the trivial cases where a — or vr/4), Axu is greater 
than AEif, rather than just greater than or equal to, as 
stated in Refs. 0, • As Axu is equal to the asymptotic 
communication capacity in a single direction Q , and the 
asymptotic bidirectional communication capacity is no 
larger than 2AEu for two qubit operations our re- 
sults show that there is an inevitable trade-off involved 
in performing bidirectional communication. That is, it 
is not possible to perform an average communication of 
Axu in both directions at the same time. 

Our results also demonstrate that larger capacities Xu 

and Axjj are obtained when the dimension of the an- 
cillas is increased above 2, in contrast to the case of the 
entanglement capacities, for which an ancilla dimension 
of 2 is sufficient. This has been shown analytically for 
En in Ref. 0|, and found numerically for AEu in Ref. 
[lOj (our numerical results also support this). 

It must be emphasized that there are limitations on 
the conclusions that can be drawn from the results due 
to their numerical nature. The capacities found here 
were found by numerical maximization, which has the 
drawback that it is possible that, in some cases, a lo- 
cal maximum may have been found, rather than the 
global maximum. However, the calculations presented 
here were performed very carefully, and in some cases 
repeated many times with different random numbers, in 
order to avoid local maxima. Therefore we are reason- 
ably confident that the results presented here are very 
close to the global maxima. 

A more serious issue is that it is not possible to per- 
form calculations for arbitrarily large ancilla dimensions. 
Although it has been proven that an ancilla dimension of 
2 is sufficient for E(j , this has not been proven for AEjj . 
Here we have found that AEjj ' is unchanged as the an- 
cilla dimension is raised from 2 to 5. This result strongly 
indicates that the values obtained for AEIj are the cor- 
rect asymptotic values. Nevertheless, it is, in principle, 
possible that larger values of AEu may be achieved for 
larger ancilla dimensions than those tested. In order to 
show conclusively that there is a strict inequality between 
Axu and AEu, it would be necessary to prove that the 
asymptotic entangling capacity is achieved for an ancilla 
dimension of 2. 

In the cases of xu and Axuj it is more difficult to es- 
timate the capacity because the capacity increases with 
ancilla dimension. For Axu the capacity only increases 
by a small amount as the ancilla dimension is increased 
from 2 to 5, indicating that the results are a good ap- 
proximation of the capacity for arbitrary ancilla dimen- 
sion. For XU2 and XU3, there are significant increases in 
the capacity with the ancilla dimension up to an ancilla 
dimension of 4, but there is no increase as the ancilla 



dimension is increased to 5. This result indicates that 
an ancilla dimension of 4 is sufficient, though further re- 
sults for higher ancilla dimensions would be required to 
be confident of this. 

In addition there is the problem that it is not possi- 
ble to perform calculations for arbitrarily large ensemble 
sizes. For the operation C/i it was found that the capac- 
ity was not increased as the ensemble size was increased 
above 2, so it is likely that an ensemble size of 2 is suffi- 
cient for this operation. For the operations U2 and U3 it 
was found that the capacities Xjj^ and Ax^^ increased 
with the ensemble size up to an ensemble size of 4. It 
is possible that the capacities may be larger for larger 
ensemble sizes; this is a topic for future research. 



APPENDIX: NUMERICAL METHODS 

The numerical techniques used to search for the max- 
ima are similar to simulated annealing 19]. For the case 

of the entanglement capacity Efj , a vector of complex 
numbers, (-0^) G C'^i'^^^^'KA+dimnB) ^ represents the ini- 
tial tensor product state. The first dimTi^ numbers are 
the coefficients representing Alice's local state, and 
the last dimTis numbers are the coefficients represent- 
ing Bob's local state, \x)b- The total state is simply 
I*) — \(I))a ® \x)b- The initial values of ipi were selected 
using a Gaussian distribution followed by normalization. 
Throughout this section, all Gaussian distributions for 
complex variables are real Gaussians multiplied by ran- 
dom phases (with a uniform distribution). 

At each step another vector of complex random num- 
bers, (Aipi), was selected via a Gaussian distribution. If 
the state represented by (V'i -I- A^i)/Af (where M is a 
normalization constant) gave a larger final entanglement 
after application of U, then the coefficients (tpi) were re- 
placed with {^i + A'ipi)/Af. This technique is equivalent 
to applying simulated annealing with a temperature of 
zero, because in no case were coefficients chosen that gave 
a smaller final entanglement. It was found that this tech- 
nique provided very rapid convergence, and using non- 
zero temperatures did not provide better convergence. 

Initially the standard deviation in the Gaussian dis- 
tribution used for the increments, a, was chosen to be 1. 
The value of a was halved each time there were 1000 con- 
secutive increments tested with none providing a larger 
final entanglement. This process was terminated when a 
fell below 10~^. At this stage the progressive changes in 
the final entanglement were on the order of 10~^^ or less. 

In the case of the capacity Xjj i a vector of complex 
random numbers, (i/'j) G c(x)dim(nA0nB) ^ ^^^^ ^^^^ ^ 
represent one initial state in the ensemble, \'4')ab- Initial 
values of 4'i were selected using a Gaussian distribution 
then normalizing. The remainder of the initial states 
in the ensemble were obtained by local unitary opera- 
tions, . These unitary operations were represented 

(k) 

by matrices of complex numbers, V^ - . The initial val- 
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ues of V^j were selected by generating complex random 
numbers with a Gaussian distribution, then using Gram- 
Schmidt orthormalization on the row vectors. The proba- 
bilities were represented by a vector of real numbers (pa;)- 
Initial values were selected using a uniform distribution 
then normalizing. 

At each step, complex random numbers for the incre- 

(k) 

ments A^i and AV^j were chosen using a Gaussian dis- 
tribution with standard deviation a, and real random 
numbers for the increments Apf^ were chosen using a uni- 
form distribution from —a/2 to (t/2. These increments 
were used to create a new ensemble with state {ipi + 
AV'i)/A/' and probabilities pk{l + Apk)/{1 +J2kPk^Pk)- 

The new operations for the new ensemble were obtained 

(k) (k) 

by adding the increments AV^j to V-)j then applying 
Gram-Schmidt orthonormalization to the row vectors. 

In this case the temperature of the simulated annealing 
was not taken to be zero. The new ensemble was selected 
if the new value for the Holveo information after appli- 
cation of U, Xncw, was larger than the previous value, 
Xoid- If the new value was lower, the new ensemble was 
selected if 

R < e(xncw-Xoid)/r^ (A 

where i? is a uniform random number between zero and 
one, and r is a tolerance equivalent to the temperature. 
The value of the tolerance r was initially taken to be 
10~^. The value of the final Holevo information was 
checked every 10000 iterations; if it had decreased, then 
the tolerance was divided by 2, and if it had increased, 
then the tolerance was multiplied by 1.1. This adjust- 
ment was found to provide reasonably rapid convergence 
to the final value. 

There were two alternative schemes used to adjust the 
standard deviation, cr, in the Gaussian distribution used 



for the increments to the ensemble. The first was similar 
to that for the entanglement, except that a was halved 
if 10000 alternatives were tested with no increase in the 
Holcvo information. The other scheme was to adjust a 
such that approximately 20% of the alternatives were ac- 
cepted. Only the first scheme was used in the initial 
part of the calculation (for the first 10^ or so increments 
tested). For later parts of the calculation both alterna- 
tives were tried. 

The numerical techniques used to calculate AE^ and 
^Xu^ were similar to those used for E^^ and Xjj'' ' though 
there are minor differences. In the case of AE^jj\ the ini- 
tial state, |V')ab was taken to be a general entangled 
state, rather than a tensor product of two local states. 
This state was represented by a vector of coefficients 

{tpi) e C»dimCHA®Wi3)^ 

The case of Axjj is rather simpler than that for Xjj ■ 
Rather than it being necessary to condider a set of uni- 
tary operations and probabilities, the ensemble was sim- 
ply represented by a set of complex coefficients for 
all of the states in the ensemble. Rather than separately 
storing probabilities, these states were not normalized, 
and the normalizations of these states were taken to be 
the probabilities. 
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